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ABSTRACT 


By  applying  the  principles  of  classical  mechanics  to  approximations  of 
aerodynamic  and  ground  interaction  forces,  a  general  three-dimensional 
debris  fragment  transport  code,  DEBRIS,  has  been  written.  Certain  field 
experiments  have  been  simulated  using  airblast  originating  in  large  H  E 
explosions  and  the  simulated  and  experimental  outcomes  have  been  compared. 
Although  simulation  of  long-distance  or  high-speed  transport  by  tumbling 
or  sliding  over  ground  surfaces  is  not  satisfactory,  short-distance 
transport,  such  as  might  be  found  following  structural  collapse,  is  suffi¬ 
ciently  realistic  to  enable  DEBRIS  and  a  previously  developed  building 
response  code,  BRACOB,  to  be  ultimately  capable  of  describing  outside  wall 
debris  patterns.  It  is  recommended  that  velocity-dependent  loss  mechanisms 
be  incorporated  in  DEBRIS  to  improve  its  simulation  of  long-distance 
transport  and  that  certain  existing  uncertainties  concerning  wall  debris 
creation  and  transport  be  lessened  by  means  of  a  full-size  or  near  full- 
size  field  experiment. 

SUMMARY 

Starting  from  initial  conditions  originating  in  wall  breakup,  the  air 
blast  drag  on  a  wall  fragment  is  calculated  from  the  simple  equation  giv¬ 
ing  the  force  in  terms  of  a  dynamic  pressure  q,  fragment  area  A  and 
empirical  drag  coefficient  Cd: 

F  =  Cd  q  A  . 

In  the  case  of  skin  fraction  drag,  q  is  based  on  the  relative  air  velocity 
parallel  to  the  surface  and  F  is  a  force  parallel  to  that  component. 

Pressure  drag  is  normal  to  the  surface  and  q  is  calculated  from  the  total 
relative  air  velocity.  In  the  regime  important  to  debris  transport  Cj  can 
be  treated  accurately  as  a  constant.  Mach  numbers  never  exceed  approximately 
0.8  and  remain  high  only  for  short  periods;  consequently  the  values  of  Cd 
appropriate  to  incompressible  flow  have  been  used. 


The  ground  is  modeled  as  (a)  a  binary  spring  compressed  by  vertical 
displacement  and  (b)  a  source  of  coulomb  friction  acting  against  horizon¬ 
tal  motion.  The  two  spring  constants  (high  during  compression,  low  during 
decompression)  have  been  varied  to  provide  coefficients  of  restitution  in 
the  range  0.0001  to  0.01;  friction  coefficients  in  the  range  0.5  to  1.1 
have  been  studied.  The  ground  force  equations  are  as  follows. 


The  spring  force  on  the  i-th  corner  can  be  written  in  terms  of  the  distance 
z^  of  the  corner  below  the  ground  as: 


F  =  k 
zi  g 


z^  |,  provided  z^<  0,  otherwise 


Fzi  =  ° 

where  the  spring  constant  k  has  one  of  two  values;  the  friction  force  in 

g 

terms  of  the  spring  force  is: 

-+  *  A  -* 

v  -  k  (k*v.) 

F  =  -  - * — — —  yF 


->  -* 
Vt  -  k  (k*Vi> 


zi 


where  V  is  the  velocity  of  the  ith  corner  relative  to  the  ground,  y  is  an 

*  A 

empirical  coefficient  of  friction  and  k  is  the  unit  vector  normal  to  the 
ground. 


Kinematics  is  based  on  two  coordinate  systems,  one  fixed  to  the  moving 
fragment  and  a  second  fixed  to  the  ground;  relative  orientations  of  these 
coordinate  systems  are  expressed  at  any  time  by  means  of  Euler  angles. 
Dynamics  of  the  center  of  mass  is  separate  from  the  dynamics  of  rotation: 
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where  R^  locates  the  center  of  mass  in  spatial  coordinates,  is  a  force 
on  the  fragment  at  the  point  located  with  respect  to  the  center  of  mass  by 
the  vector  r^  and  u)  is  the  angular  velocity  about  the  center  of  mass.  Moment 
of  inertia  tensors  considered  so  far  have  been  diagonal  but  not  necessarily 
otherwise  symmetrical. 


In  integrating  equation  (1)  a  centered  scheme  is  used: 

?  d2MO 

R  (t  .  .)  =  At  - +  2R  (t  )  -  R  (t  .) 

c  '  n  +  1  ^2  c  n  c  n-1 

but  because  of  discontinuities  in  the  Euler  angles,  equation  (2)  for  orienta¬ 
tion  is  solved  by  an  uncentered  integration  scheme: 


AS$  (t  )  =  w(t  )  At  +  w(t  )  At^ 


where  At  is  the  time  increment  and  A^  is  the  angular  displacement  of  the 
fragment  during  the  increment. 


The  realism  of  the  three-dimensional  model  (the  FORTRAN  code  of  which  is  named 
DEBRIS)  has  been  tested  against  experimental  observations  at  PRAIRIE  FLAT  and 
DIAL  PACK.  At  low  overpressures  and  for  transport  over  short  distances,  the 
model  xs  very  good.  Simulated  individual  concrete  masonry  units  (CMU)  blown 
out  of  a  wall  at  PRAIRIE  FLAT  by  a  9  psi  wave  appear  to  travel  very  nearly 
the  same  distance  as  did  the  actual  masonry  units;  the  simulation  of  a 
tumbling  plywood/concrete  cube  at  DIAL  PACK  can  be  accurately  simulated  when 
peak  overpressure  is  15  psi.  By  adjustment  of  the  coefficient,  the  simulation 
of  the  tumbling  cube  at  30  psi  can  be  made  satisfactory.  However,  at  higher 
pressures  (e.g.,  50  and  100  psi)  the  simulation  clearly  shows  the  need  of 
additional  loss  mechanisms  than  coulomb  friction;  a  viscous  or  velocity- 
dependent  force  has  been  suggested  to  overcome  this  deficiency. 


Another  kind  of  experiment  has  been  simulated  with  both  the  transport  model 
and  the  previously  developed  code  BRACOB  (Blast  Response  and  Collapse  of 
Buildings).  At  PRAIRIE  FLAT  a  small  masonry  building,  10  ft.  on  a  side, 
was  exposed  to  two  peak  blast  pressures,  9  and  30  psi,  and  the  final  disposition 


3 


of  wall  debris  was  recorded  photographically.  Simulating  the  blast  impact, 

BRACOB  predicts  a  sequence  of  wall  displacements  and  failures.  With  initial 
or  departure  conditions  chosen  by  means  of  BRACOB  and  by  test  results  clearly 
showing  breakup  patterns,  the  transport  model  can  predict  the  impact  points 
and  final  locations  of  wall  fragments.  The  predicted  scenario  at  30  psi 
appears  relatively  unambiguous,  and  observed  wall  debris  patterns  are  consis¬ 
tent  with  predicted  patterns;  apparently  edge  support  gives  rise  to  arching 
or  bowing  of  the  panel,  which  in  turn  determines  the  direction  of  departure. 
Fragment  interaction  takes  place  just  downwind  of  the  building. 

At  9  psi  the  observed  pattern  is  not  symmetrical,  probably  due  to  slight 
assymmetries  in  blast  or  construction  and,  moreover,  the  simulation  is  not 
as  clear  as  at  30  psi  because  of  the  predicted  absence  of  classical  air- 
blast  flows  through  the  building.  However,  simulation  and  observation  are  not  in¬ 
consistent.  The  two  codes  BRACOB  and  DEBRIS  appear  to  promise  means  of 
relating  wall  debris  trajectories  and  final  patterns  to  events  during  wall 
collapse. 

In  preparation  for  making  estimates  of  urban  debris,  the  simulation  system 
BRACOB  has  been  greatly  enhanced.  It  can  now  handle  reinforced  concrete 
wall  panels  attached  to  yielding  frames  as  well  as  reinforced  walls,  and 
storage  has  been  increased  to  make  analysis  of  office  buildings  and  hospitals 
possible.  A  test  using  construction  drawings  of  a  Philadelphia  hospital  shows 
the  enhanced  code  fully  operational. 
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I  INTRODUCTION 


Among  the  effects  of  nuclear  weapons  upon  an  urban  area,  those  produced 
by  fire  are  probably  the  least  understood,  yet  fire  is  undoubtedly  one 
of  the  most  awesome  effects.  Fire  researchers  interested  in  nuclear  explosions 
encounter  many  difficulties  in  trying  to  understand  fire  ignition  growth  and 
spread,  not  the  least  of  which  are  the  difficulties  stemming  from  the  effects 
of  blast:  (1)  the  direct  effect  of  blast  wind  and  pressure  on  fire  and 
(2)  the  indirect  effects  of  blast  through  its  creation  of  debris.  (Debris 
is  here  defined  to  be  any  solid  object  displaced  by  blast  from  its  normal 
position  or  function.)  It  is  the  second  of  these  two  aspects  of  the  fire- 
blast  interaction  that  inspired  this  research. 


Thermal  irradiation  of  the  urban  complex  begins  at  the  moment  of  weapon 
detonation.  After  only  an  inconsequential  delay,  the  blast  wave  also 
departs  the  fireball.  Because  its  travel  to  the  irradiated  objects  is  much 
slower  than  that  of  the  radiation,  at  large  distances  from  the  explosion  the 
two  effects  are  readily  separable  in  time,  but  at  small  distances  the 
temporal  separation  is  not  as  sharp.  This  statement  is  made  quantitative  in 
Figure  1,  where  the  contours  of  thermal  fluence  are  shown  as  dashed  lines 
in  coordiantes  of  time  on  the  ordinate  and  range  or  overpressure  on  the 
abscissa.  Since  the  thermal  pulse  shape  is  the  same  at  all  ranges,  the  time 
axis  is  also  marked  to  show  certain  critical  times  during  irradiation,  i.e., 
the  times  when  fluence  is  20%,  50%  and  80%  of  the  total.  The  heavy  solid 
line  starting  at  the  lower  left  and  curving  to  the  upper  right  corner  repre¬ 
sents  the  time  of  blast  arrivals  as  a  function  of  range.  Two  dissimilar 
regions  of  overpressure  are  highlighted  with  shading.  At  ranges  correspond¬ 
ing  to  peak  free-field  overpressures  near  50  psi  (345. kPa),  only  20%  of  the 
thermal  fluence  has  arrived  by  the  time  the  blast  begins  to  create  debris; 
whereas  at  overpressures  near  4  psi  (27.6  kPa),  blast  arrival  follows  the 
arrival  of  80%  of  the  thermal  fluence.  Furthermore,  assuming  20  cal/cm^  is 
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the  threshold  for  primary  ignitions  at  the  4  psi-(27.6  kPa)  range,  such 

1/3 

ignitions  will  have  well  over  5  s/mt  to  develop  before  experiencing 
blast  effects;  comparable  development  time  at  50  psi  (345  kPa)  on  the  other 
hand,  is  less  by  more  than  a  factor  of  ten. 

As  time  after  blast  arrival  passes  at  a  point  in  an  urban  complex,  the  pro¬ 
blem  of  accurately  tracing  debris  transport  and  its  interaction  with  other 
debris  and  other  structures  becomes  almost  impossibly  complicated.  However, 
immediately  after  blast  impact,  the  problem  is  relatively  tractable  and  the 
present  study  begins  at  that  point.  After  certain  necessary  theoretical  and 
experimental  preliminaries,  this  report  addresses  the  problem  of  the  conver¬ 
sion  of  the  exterior  walls  of  a  building  into  debris  and  its  subsequent 
transport. 

The  physical  conditions  determining  collapse  of  the  exterior  walls  of  a 
building  under  airblast  loading  have  been  described  elsewhere  (1,  2).*  The 
physics  of  transport  of  such  debris  are  discussed  in  the  next  section. 


*  References  are  included  at  the  end  of  the  report. 
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II  PHYSICAL  PRINCIPLES  OF  DEBRIS  TRANSPORT 


Forces  transporting  debris  fragments  come  from  two  sources:  (1)  blast  wind 
and  (2)  impact  with  the  ground,  walls,  or  other  objects.  The  following  dis¬ 
cussion  of  wind  forces  and  ground  forces  presents  a  calculation  method  of 
predicting  debris  transport. 

Wind  Forces 

Generally,  the  translating  effects  of  pressure,  as  distinct  from  the  wind, 
are  negligible  on  a  debris  fragment  because  the  so-called  "static"  pressure 
is  short-lived  against  all  but  large  objects  such  as  building  facades;  debris 
transport  for  the  most  part  is  effected  by  the  "dynamic"  pressure,  q,  defined  as 

q  =  h  pu2  (1) 

where  p  =  local  air  density 
and  u  =  air  particle  speed. 

Forces  on  ohjects  may  he  expressed  in  nondimensional  coefficient 
form  in  terms  of  lift,  drag,  and  side-free  coefficients.  The  force  in 
the  windward  direction  is  given  by  the  drag  coefficient,  the  dynamic 
pressure,  and  a  characteristic  reference  area  for  the  body,  as 

F  =  CdqA  .  (2) 

In  this  expression  F  is  the  force  in  the  windward  direction,  Cd  is  the  drag 
coefficient,  and  A  represents  the  characteristic  area  of  the  body.  Often  A  is 
taken  as  the  cross-sectional  area  of  the  object  projected  normal  to  the  upst.eam 
flow  direction.  The  drag  may  be  divided  into  components.  "Pressure"  drag  is 
created  by  the  local  fluid  forces  normally  exerted  at  the  surface.  "Friction" 
drag  is  created  by  tangential  fluid  friction  at  the  surface.  In  magnitude,  the 
friction  drag  is  generally  of  the  order  1/10  of  the  pressure  drag.  Obviously, 
the  local  normal  pressure  forces  on  a  surface  can  create  both  drag  and  lift 
components,  depending  on  the  orientation  of  the  surface  to  the  upstream  wind 
vector,  as  illustrated  in  Figure  2. 
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Let  us  consider  a  body  made  up  of  flat  surfaces.  The  body  may  be  divided 
into  two  parts:  (1)  the  part  that  faces  toward  the  air  stream  (fore  body) 
and  (2)  the  part  that  faces  downstream  (base  region).  The  drag  coefficient  may 
be  divided  into  fore  body  drag  and  what  is  commonly  called  base  drag  to 
account  for  the  separate  influences  of  the  two  parts  of  the  body.  In  the 
present  work,  base  drag  and  fore  body  drag  have  been  combined.  Although 
the  pressure  on  a  flat  surface  on  which  the  flow  is  impinging  is  usually 
almost  uniform  over  the  surface,  the  pressure  over  a  trailing  flat  edge  may 
not  be  uniform.  The  resulting  forces  can  contribute  a  net  moment  to  the  body 
while  it  is  in  flight,  depending  upon  its  orientation  to  the  flow.  The 
moment  may  result  in  stable  or  unstable  flight,  depending  upon  the  location 
of  the  body's  center  of  gravity. 

Numerical  values  used  to  estimate  the  drag  coefficient,  C^,  were  taken  from 
wind  tunnel  data  that  covered  a  range  of  shapes  and  wind  incidence  angles 
(  3,  4  ).  The  range  of  shapes  and  angles  actually  tested  and  reported  in 
the  literature  is  limited,  but  there  are  enough  similarities  between  tested 
bodies  and  debris  fragment  shapes  to  make  extrapolation  reasonable.  Also, 
computational  simulations  of  reported  observations  of  actual  building  debris 
transport  serve  as  a  rough  check  of  the  drag  coefficients  used. 


Both  friction  and  pressure  drag  coefficients  depend  on  the  dimensionless 
Reynolds  number. 


R 

e 


V  p  L 
r  y 

1 


where 


and 


relative  velocity  of  air  and  object 
characteristic  dimension  of  object 
physical  vicosity  of  air. 


Figures  3A  and  3B  illustrate  this  dependence  of  drag  on  Reynolds  number. 
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Circular  and  Square  Plates  in 
Normal  Flow 


The  pressure  drag  in  (A)  of  Figure  3  is  for  a  flat  two-dimensional  or 

square  plate  normal  to  the  flow.  The  pressure  drag  curve  is  typical  in 

that  it  shows  high  drag  at  low  Reynolds  numbers,  but  becomes  constant 

at  high  numbers.  The  Reynolds  numbers  for  the  flows  considered  in  this 

4  6 

report  are  generally  in  the  range  10  to  10  .  Therefore,  the  pressure 
drag  coefficients  will  be  nearly  constant  with  Reynolds  number. 

Drag  coefficients  have  been  measured  in  wind  tunnels  for  a  number  of 
regular  shapes  such  as  cubes,  spheres,  plates  and  cylinders.  For  example, 
such  pressure  drag  information  exists  for  a  flat  plate  as  a  function 
of  flow  incidence  angle  as  shown  in  Fig,  4.  Cuhes  placed  on  a  wind 
tunnel  floor  have  a  drag  coefficient  of  1.05  with  one  face  normal  to 
the  flow  and  0.90  when  tilted  at  45°.  The  lower  value  for  a  cube  is 
approximately  twice  the  coefficient  for  a  sphere  (0.47)  (3). 

Drag  coefficients  have  also  been  deduced  for  some  of  these  same  shapes 

during  drop  tests  in  the  atmosphere  from  the  time  of  fall,  with  the 

assumption  the  coefficient  was  constant  during  the  fall  (5) .  For 

example,  cubes  dropped  with  one  face  down  yielded  values  of  in  the 

range  1.07  to  1.3.  Terminal  Reynolds  numbers  for  these  measurements 

4  4 

were  in  the  range  2  x  10  to  8  x  10  .  The  study  was  extended  to  include 
"irregular"  plates  of  random  shape  (i.e.,  fractures  glass  plates); 
statistically,  was  found  to  differ  from  that  for  "regular"  plates  by 
no  more  than  5.6%,  which  was  within  the  scatter  of  data. 

Using  various  three-dimensional  shapes  made  from  lightweight  wood.  Bell 
Telephone  Laboratories  measured  average  drag  coefficients  by  means  of 
what  were  essentially  drop  tests  in  a  vertical  wind  tunnel  (6). 

The  largest  dimension  of  the  objects  varied  from  8  to  24  inches  (203  to 
to  610  mm),  and  the  least  wind  velocity  needed  to  hold  the  object  in 
suspension  was  recorded.  Although  measured  drag  coefficients  varied 
widely,  they  were  all  of  the  order  of  magnitude  of  unity.  For  example, 
a  balsa  wood  cube,  8  inches  (203  ram)  on  edge,  showed  an  average  drag  coeffi¬ 
cient  of  1.5,  as  did  a  simulated  cinder  block  made  of  the  same  material. 


8 


The  drag  coefficients  used  in  the  present  study  were  taken  from  wind- 
tunnel  test  data  for  "comparable"  shapes,  because  good  agreement  existed 
between  the  results  of  drop  tests  and  wind  tunnel  measurements  of  drag 
coefficients.  For  example,  masonry  units  were  treated  as  cubes;  large 
wall  fragments  were  considered  to  be  flat  plates.  The  drag  coefficients  of 
a  tumbling  cube  were  approximated  to  be  the  average  between  the  value  with 
the  leading  face  normal  to  the  flow  (1.05)  and  the  value  with  two  leading 
faces  45%  to  the  flow  (0.80).  The  value  of  reference  area  A  used  in 
equation  (2)  is  the  maximum  cross-sectional  area  of  the  body  in  a  plane 
normal  to  the  relative  wind  velocity.  The  drag  coefficients  for  flat  plates 
are  given  in  Figure  4.  They  were  corrected  for  Reynolds  number  using  the 
values  from  (a)  of  Figure  3.  Friction  drag  was  included  for  all  leading 
surfaces  using  the  values  from  (B)  of  Figure  3.  The  air  viscosity  was 
computed  based  upon  the  air  temperature  at  the  fragment  surface. 

Ground  Forces 

In  an  attempt  to  model  torques  arising  from  ground  impact,  all  calculated 
ground  forces  were  considered  applied  to  corners  of  the  debris  fragment. 

When  by  symmetry  torques  arising  from  different  corners  should  cancel, 
special  numerical  techniques  were  employed  to  ensure  they  did  so.  Energy 
dissipation  was  provided  by  representing  the  ground  as  a  single  composite 
spring,  relatively  stiff  in  compression  and  weak  in  decompression  (7). 

Values  chosen  for  spring  constants  corresponded  to  a  coefficient  of  resti¬ 
tution  equal  to  .01.  Spring  compression  began  at  the  moment  the  corner 
impacted  the  ground,  after  which  the  spring  applied  increasingly  upward- 

directed  force  on  the  corner  until  its  downward  motion  reversed.  At  that 

-4 

time,  the  spring  constant  was  reduced  by  a  factor  of  10 

A  second  ground  interaction  is  through  sliding  friction,  here  modeled  as 
coulomb  friction,  proportional  in  magnitude  to  the  force  in  the  ground 
spring  and  a  friction  coefficient  y.  Its  direction  is  opposite  to  the 


component  of  fragment  velocity  parallel  to  the  ground  surface,  i.e.. 


V  +  V 
rx  ry 


x  i  +  y  j 


V  +  V 
rx 


nl 


ufz=  - 


P  +  ?! 


. ,  l/2  7 


This  expression  is  written  in  a  coordinate  system  in  which  the  Z-axis 
is  normal  to  the  ground  surface.  Values  of  y  examined  in  this  work 
varied  from  0.50  to  1.1. 
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Ill  EQUATIONS  OF  MOTION 


Because  debris  transport  is  not  simple  downwind  translation,  three- 
dimensional  equations  of  motion  have  been  written  to  take  into  account 
initial  displacement  and  rotation,  ejection  into  crosswinds,  and  ground 
impact  forces  and  torques.  Not  yet  accounted  for  are:  fragment  breakup 
in  flight  or  on  impact,  interactions  between  fragments,  non-ideal  airflows 
such  as  those  undoubtedly  existing  at  least  transiently  around  a  collapsing 
building,  and  eccentric  wind  forces  on  flat  plates. 

Coordinate  Systems 

Throughout  the  following  discussion,  reference  is  made  to  two  coordinate 
systems:  the  first  fixed  in  space,  with  the  Z-axis  normal  to  the  ground 
surface  and  the  Y-axis  normal  to  the  wall  in  its  undisturbed  position  or 
to  one  face  of  the  fragment,  and  the  second  fixed  to  the  fragment  or  cbject 

*  n  A 

transported  in  which  unit  vectors  are  a,  p,  and  y. 

A  A 

Before  the  fragment  is  disturbed,  Y  coincides  with  the  unit  vector  k  along 
/\  ^  A’ 
the  Z-axis,  a  with  the  unit  vector  i  along  the  X-axis,  and  B  with  j,  the 

unit  vector  along  the  Y-axis.  Solution  of  the  equations  of  motion  will 

then  provide:  (1)  the  instantaneous  relative  orientation  of  the  two 

coordinate  systems  and  (2)  the  position  of  the  center  of  fragment  mass  as 

a  function  of  time. 

Any  angular  displacement  of  a  rigid  body  can  be  expressed  as  a  sequence 
of  the  displacements  through  the  Eurler  angles  0,  0,  and  t p,  illustrated  in 
Figure  5.  (8).  Once  the  Euler  angles  are  known,  any  vector  fixed  in  the 
fragment  can  be  expressed  in  the  spatial  coordinate  system  by  a  simple 
tensor  transformation: 


^•1 
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T  =  M  Mr.  M , 

Y  o  4> 


where 


'  cost#  sin<#  0 

M  =  / 

<#  l-sin<#  cost#  0 


(3a) 


0  0 
cos0  sin0 
-sin©  cos0 / 


(3b) 


M 


/cost#  sini#  0 


^  |-sin<#  cost#  0  |  (3c) 

0  0  1 

Tensor  T  contains  the  components  of  the  unit  vectors  explicity: 


T  = 


ka 

jS  kB 

\  V 


°k 


Dynamics. 

The  physics  of  the  calculation  is  contained  in  the  following  equations. 

Motion  of  the  center  of  mass  follows  from: 


■  -  m 


w 


space 


where  11  is  the  position  of  the  center  of  mass  in  the  spatial  coordinates. 


M  =  mass  of  body 

F  =  E  the  sum  of  all  forces  on  the  body 
t  =  time. 

The  subscript  "space"  means  R  is  expressed  in  the  spatial  frame. 

The  following  equations  give  angular  displacement  as  a  function  of  time: 
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*  =(£) 


space 


+  li)  *  L 


where  L  =  angular  momentum  about  center  of  mass  (CM) 

0)  =  angular  velocity  about  CM 
N  =  I  sum  of  all  torques  about  CM  , 

Angular  momentum  of  a  rigid  body  is  the  tensor  product  of  angular  velocity 
and  a  moment  of  inertia  tensor  1 : 


L  =  I  W 


Moment  of  inertia  is  a  property  of  the  distribution  of  mass  in  the  body,  (see 
reference  (3)  for  details  of  calculation)  but  for  all  cases  considered  so  far, 
X  is  diagonal,  i.e., 

10  0 
1=0! 
o  rB  o 

0  0  L. 


in  which  the  values  of  the  components  are  here  treated  as  constants.  The 
derivative  of  angular  momentum  becomes: 

($)**-  +  +  <« 

and  u  x  L-oi  *  I'd  -  (Iy  -  IgH  ^{1^  -  1^)  (7) 


+  tve'iff  S' 


where  u)  ,  a)R,and  w  are  time  derivatives  of  the  components  of  angular 
velocity  in  the  body  coordinate  system.  Addition  of  the  tern  to  x  L  accounts 
for  the  change  of  unit  vectors  a,  S,  and  y  with  time. 


The  torque  N  can  be  written: 
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N  =  Z  r±  x  t±  =  i  {  a  (rigFiY  -  r^HBO:^  F^-rJ^) 


+  y(r .  F. D  -  r . „F.  )  } 
'  ia  iB  iB  ia 


(8) 


where  is  the  vector  from  the  CM  to  the  point  of  application  of  the  i-th 
external  force  F^,  and  the  summation  is  over  all  external  forces.  Dropping 
the  subscript  i,  any  force  on  the  fragment  can  be  expressed  in  body 
coordinates  through  the  relations  : 


F  =  Fa  +  F  a  +  F  a 
a  xx  y  y  z  z 

FQ  =  F  8  +  F  8  +  F  B 

B  xx  yy  z  z 

F  =  F  y  +  F  B  +  F  B 
Y  xx  y  y  z  z 


Integration 

Since  the  external  forces  are  readily  calculable  in  terms  of  the  spatial 
coordinate  system,  the  equations  of  motion  of  the  CM  are  integrated  directly 
from  the  approximation  to  the  second  derivative,  e.g.: 

R  (t  -  R  (t  )  R  (t  )  -  R  (t  -1) 
x  n+1  xN  n _ x  n  x  n  / 

At  At 


dt2  At 


or 


d2  R 


R  (t  ..) 
x  n+1 


=  (At) 


x 


dt 


,  t  2R  (t  )  -  R  (t  .) 
2  x  n  x  n-1 


(9) 


In  these  expressions  Rx  is  the  x-component  of  the  location  vector  R  of  the 
CM  in  the  spatial  coordinates;  the  argument  t  identifies  the  value  of  the 
displacement  at  the  time  of  the  n-th  time  step;  and  At  is  the  magnitude  of 
the  time  step,  which  is  constant.  Equation  (9)  provides  a  value  of  one 
component  of  CM  displacement  in  terms  of  previously  calculated  values.  There 
are  analogous  equations  for  the  other  two  components. 
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The  value  of  the  second  time  derivative  in  equation  (4)  follows  from 
equation  (3) . 

d  Rx  1  i  F  . 

— =  —  xi 

dt  M 

The  unit  vectors  in  the  spatial  coordinate  system  are  constant  in  time. 

Calculation  of  the  fragment  orientation  consists  of  two  steps: 

(1)  Determination  of  the  infinitesimal  rotation 

AS^  =  to  At  +  1*  to  (At)^  (10) 

and  (2)  computation  of  the  Euler  angles  expressing  the  new  orientation 
at  tR  +  Equations  of  motion  can  be  written  directly  in  terms  of  first 
and  second  derivatives  of  the  Euler  angles,  but  the  angles  and  their  first 
derivatives  are, in  general, strongly  discontinuous,  requiring  repeated 
changes  of  coordinate  systems  during  numerical  integration.  Therefore, 
although  the  scheme  represented  by  equation  (10)  employs  uncentered  time 
steps,  it  is  far  less  cumbersome  than  an  analogue  of  equation  (9).  During 
actual  numerical  integration, the  accuracy  of  equation  (10)  was  tested  by 
examining  the  results  for  conservation  of  angular  momentum. 

The  derivatives  in  equation  (10)  are  found  from  equations  (5),  (6), and 

(7)  after  the  instantaneous  value  of  the  torque  N  is  calculated  in  equation 

(8)  i.e., 


IN  -  0)DW  (I  -  I0)l  /i 

La  3  Y  Y  &  a 

N0  -  (l)  U>  (I  -  I  )1  / Ia 

3  Y  a  a  y  J  '  3 

N  -  u)  u)„  (I„  -  I  )]  /l 

Y  a  3  3  a  '  y 


VW  ’  W  +  4t“s(tn) 

YW  ’  VV  + 


The  displacement  A£l  is  an  infinitesimal  rotation  through  the  angle  | Ai2 1  about 
an  axis  through  the  CM  parallel  to  A£^.  Such  a  displacement  can  be  expressed 
by  Euler  angles  <}>',  0^,  and  tp'  in  the  following  way  (see  Figure  6).  The  first 

|  A  A 

angle  <p  is  a  rotation  of  the  body  coordinate  system  about  y  to  align  g 
with  the  projection  of  AS^  in  the  ag  plane.  The  second  angle  8  (see  Figure  5) 

/s  ^  -j*- 

is  a  rotation  of  coordinates  about  the  displaced  axis  a  aligning  y  with  AQ, 

/  A 

and  finally  4*  ,  the  Euler  rotation  of  the  body  about  the  displaced  axis  a,  is 

chosen  so  that 

ip'  =  |  VQ  ( 

/V  /V  A 

The  coordinate  axes  a,  g,  and  y  are  returned  to  their  original  positions  in 

A  /v 

space  by  inverse  rotations,  first,  -  0  about  the  displaced  a  and,  second, 

— <p  about  displaced  y.  The  transformation  is  expressed  in  terms  of  tensors 
similar  to  those  defined  in  equation  (3) : 


where 


"e 


v  vv 


-1  *  2  u 

Q'  =  tan  (Afl  a  +  Aflp)^ 
A(2y 

|  A  i?  | 


when 


An  =  A no  =  0,  then 
a  8  ’  r 


:  0^=  0. 


The  final  step  is  the  discovery  of  the  Euler  angles  0,  0  ,  and  ij;  describing 
the  new  orientation  as  a  displacement  from  the  undisturbed  position,  and 
this  is  done  through  the  tensor  equation: 

Wr  TW  (tij)  (n> 

where  <J>,  0  and  ^  are  the  Euler  angles  at  the  preceding  time  step.  Since 
is  known,  equation  (11)  leads  to  a  series  of  equations  that  can  be 
solved  in  the  following  order: 
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COS  0  = 

cos  <}>  =-t32^  sin  ® 
sin  <j>  =-t^/sin  0 
cos  tp  =  C^^/sin  0 
sin  iji  =  t^/sin  Q 


when  sin  0=0  then 


cos  (  +  ip  ) 


C11  +  t22 


The  values  <}> ,  0,  and  ip  become  4>>  0  and  ^  for  the  next  integration  step. 

The  foregoing  calculational  scheme  has  been  embodied  in  the  FORTRAN  code 
DEBRIS 

Incidentally,  the  values  of  4>,  9,  and  $  calculated  as  above  could  be  used  to 
identify  discontinuities  in  the  Euler  angles  in  the  centered  integration  scheme 
outlined  below. 

The  method  is  applied  to  the  second  time  derivatives  of  the  Euler  angles, 

(j),  0,  and  which  arise  from  the  derivative  of  the  angular  velocity  ta,  as 

■~y 

follows.  The  infinitesimal  displacement  <a  At  may  be  represented  by  the  vector 
sum  of  the  three  Euler  rotations: 

A  A-  All 

a)  At  =  k  A$  +  i  A9  +  k  Aiji 
or,  as  long  as  A<p,  A0  and  A<p  are  infinitesimals, 

-*■  .A>  .A" 

u>  =  0k  +  0i  +  4»k  .  (12) 


But  in  terms  of  the  body  coordinate  system: 


=  SIN  ip  SIN0  a  +  COS  ^  SIN  9  &  +  COS  9  y 


COS  ^  a  -  SIN  ^  3 
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so  that 


u>a  »  <p  sin  ip  sin  0  +  e  cos  ip 

0Jo  =  <p  COS  ip  SIN  0  -  0  SIN  ip 

p 

oj^  =  |  cos  e  +  ip 

and 

•  ••  •• 

«*)  =<>  SIN  ip  SIN  0  +0  COS  ip 

+  0  ip  COS  ip  SIN  0  +  |  0  SIN  ip 

-  0  ip  SIN  ip 

_  •«  »• 

Wg  =  <p  COS  Ip  SIN  0  -  0  SIN  ip 

-  |  Ip  SIN  Ip  SIN  0  +  <p  0  COS  tp 

-  0  ip  COS  Ip 

«•  •• 

U)^  =  <p  COS  0  +  ip  -  <p  0  SIN  0  . 

The  three  equations  resulting  from  substituting  these 
equation  (5)  are: 

<p  I  SINlp  SIN0  +  01  COSip 
ct  ct 

-  N  -  I  {  $  §  COS0  SINlp  -  Ip  4 
+  0  ip  SINlp}  +  (Ig  -  I y){yz4> 

-  ip  0  SINlp  +  <p  ip  COSip  SINS  - 


COS  0 

COS  0 


expressions  in 

'  COSip  SINS 
COSip  SIN20 
<p  0  COS0  SINlp} 
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<f>  COSip  SIN0  -  0  Ig  SINip 

=  N0  -  I0{-<f  if  SINlp  SINS  +  <f  0  COS0  COSip  -  0  if  COSip} 
P  P 

+  (I  -  I  )  {  ‘/2  i2  SIN20  SINip  +  <f  0  COSd  COSip 
Y  ci 

+  if  f  SINip  SIN0  +  if  0  COStp} 


d>  I  COS0  +  ip  I 
Y  Y 

=  N  +  i  je  SIN0  +  (I  -  Ift)  {  Vz  f2  SIN2ip  SIN20 

IT  (X  p 

.  .  , 

+  ij>  0  C0S2ip  SIN0  -y20  SIN2lp}  . 


These  equations  have  solutions  for  <J>,  0,  and  ip,  as  long  as  SIN0  =£■  0. 

For  that  reason,  coordinate  systems  (which  coincide  initially)  may  be 
chosen  such  that  the  initial  orientation  of  the  fragment  corresponds  to 
0^-0,  tt,  2v„  However,  during  the  subsequent  calculation, discontinuities 
in  the  Euler  angles  may  arise,  as  can  be  understood  by  considering  the 
sudden  onset  of  a  torque  about  the  3-axis.  The  Euler  system  as  presently 
constituted  can  account  for  an  infinitesimal  rotation  about  3  only 

by  first  introducing  an  exchange  of  coordinates  (i.e.,  increase  ^  or  0 

u 

discontinuously  by  i-?).  In  other  words  A<j>,  A6,  and  Alp  in  equation  (12) 
must  be  kept  infinitesmal  by  redefining  the  Euler  angles  <J>,  0,  and  ip., 

Test  Cases 

To  test  the  realism  of  the  simulation,  a  number  of  trajectories  were 
calculated.  Every  case  was  paralleled  by  an  actual  incident  of  full-scale 
debris  transport  under  loading  by  a  blast  stemming  from  an  HE  explosion 
equivalent  to  1  kt  (4.2  TJ)  nuclear  (HOB  =  0),  for  which  at  least  partial 
data  exists. 


Single  Concrete  Masonry  Unit  (CMU) 


Two  CMU  walls  (one  of  which  is  sketched  in  Figure  7)  were  exposed  head- 
on  a  9-psi  blast  at  operation  PRAIRIE  FLAT  (6).  Although  one  walj.  was 
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mortared  and  one  was  not,  the  debris  patterns  were  very  similar.  Each  was 
5  1/3  ft  (1.62  m)  high  and  12  ft  (3.65  m)  long.  Using  the  mass  distribution 
in  a  typical  U. S.  two-cell  CMU  weighing  27  lbs  (12.3  kg)*,  and  assuming 
zero  initial  velocity,  a  unit  in  the  top  course  is  transported  by  blast  wind 
a  maximum  distance  of  13.9  ft  (4.24  m)  as  shown  in  Fig.  7  ,  when  p  =  0.5. 
Changing  p  to  1.0  results  in  a  trajectory  ending  at  18.5  ft  (5.64  m) .  The 
greater  travel  results  from  the  greater  torque  given  the  unit  while  in 
sliding  contact  with  the  ground.  Apparently  the  relatively  large  torque 
has  the  effect  of  increasing  the  height  of  the  bounce  and  hence  the  length 
of  travel.  As  will  be  seen  later,  increasing  the  coefficient  of  friction  p 
increases  energy  dissipation  until  p  becomes  unity,  at  which  point  total 
dissipation  starts  to  decrease.  Although  no  measurements  of  final  debris 
locations  are  reported,  the  post-shot  photograph  in  Figure  9  shows  maximum 
downwind  displacements  consistent  with  that  calculated  for  p  -  0.50. 


*  Moment  of  inertia  in  CGS  units  is: 

0.95535  x  106  0. 

0.  0.21975  x  107 

0.  0. 


0. 

0. 

0.23326  x  10' 
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Figure  7  Preshot  View  of  CMU  Wall  ,  PRAIRIE  FLAT 
Pso  =  9  psi  (62  kPa). 


I 
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Trajectory  of  center  of  mass 


Horizontal  Distance  (ft) 


Figure  9  CMU  Wall  Postshot  View  of  PRAIRIE  FLAT 

P  9  psi  (  62  kPa  ) 


In  another  PRAIRIE  FLAT  experiment,  a  four-walled  building  approximating 
a  cube  10  ft  (3.05  m)  on  edge  (Figure  10)  was  exposed  head-on  at  two 
overpressures,  9  and  30  psi  (62.0  and  207  kPa) .  Although  the  collapse  of 
these  structures  will  be  discussed  later,  the  simulated  flight  of  a 
single  CMU  from  the  sidewall  will  be  used  here  to  demonstrate  conservation 
of  angular  momentum  and  indicate  the  power  of  debris  analysis.  A  unit 
in  the  top  course  was  assumed  to  depart  the  wall  54  ms  after  a  30-psi 
(207~kPa)  blast  impact  with  a  velocity  of  magnitude  69  ft/s  (21.0  m/s) 
directed  normally  to  the  blast  wind.  The  calculated  trajectory  of  the  CM 
is  represented  in  Figures  11A  and  11B.  Figure  11A  presents  the  projec¬ 
tion  of  the  trajectory  in  the  plane  normal  to  the  ground  and  parallel  to 
the  wind  while  Figure  11B  gives  the  projection  of  the  trajectory  in  the  plane 
normal  to  the  wind.  The  CMU  was  still  moving  after  a  crosswind  travel  of 
S5  it  (25.9  m)  and  travel  of  24  ft  (7.315  m)  dowuwind .  Wind  reversal 
occurs  before  the  unit  strikes  the  ground  but  does  not  reverse  the  down¬ 
wind  CMU  motion  until  near  the  end  of  the  flight.  Ground  friction  p  in 
this  case  was  0.90. 

Conservation  of  angular  momentum  was  tested  after  the  second  ground  impact 
at  1.4301  s  or  6,621  calculational  cycles  after  the  first  ground  impact  at 
0.7680.  Results  are  tabulated  in  Table  1. 

Table  1 

CONSERVATION  OF  ANGULAR  MOMENTUM 


Time  (s) 

W  (cgs) 

Ly 

■ _ lat 

1.4301 

-  — £ 

-0.5138  x  10 

0.3710  x  107 

-0.4081  x  10 

1.5301 

-0.5144 

0.3666 

-0.4106 

1.6301 

-0.5150 

0.3692 

-0.4057 

1.7301 

-0.5156 

0.3723 

-0.4113 

1.8301 

-0.5163 

0.3670 

-0.4104 

1.9301 

-0.5169 

0.3723 

-0.4076 

Constancy  of  angular  momentum  is  probably  good  enough  for  present 
purposes  and  a  major  effort  to  improve  the  integration  accuracy  is  not 
justified. 
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mulated  Trajectory  of  Single  CMU  from  Top  of  Sidewall  ( Downwind) 


An  interesting  plot  of  the  calculated  angular  displacement  of  the  vector 
^  during  the  interval  1.78  to  1.92  s  is  shown  in  Figure  12.  The  projection 
of  the  vector  in  the  x  -  y  plane  is  marked  by  the  dashed  line.  The 
lack  of  symmetry  in  the  inertia  tensor  prevents  true  precession. 

Starting  the  CMU  with  the  same  initial  speed  and  departure  time,  but  half¬ 
way  up  the  sidewall  instead  of  in  the  top  course,  and  allowing  it  to  come 
to  rest  before  stopping  the  calculation,  results  in  a  final  downwind  dis¬ 
placement  of  19  ft  (5.79  m)  and  crosswind  displacement  of  93  ft  (28.3  m) 
after  4.47  s  of  travel  (Figures  13A  and  B).  As  the  postshot  photograph 
in  (figure  14B)  shows,  this  appears  to  be  an  unrealistically  large  amount 
of  travel,  suggesting  either  the  sidewall  does  not  break  up  by  ejecting 
CMU  directly  into  the  crosswind  or  that  the  initial  conditions  are  in  error. 
For  example,  changing  the  initial  conditions  to  a  start  at  54  ms  with  a 
speed  of  34  ft/s  (10.4  m/s)  results  in  net  displacements  of  11.6  ft  (3.53  m) 
downwind  and  25.6  ft  (7.80  m)  crosswind,  which  are  clearly  more  realistic. 

In  both  these  examples,  the  unit  had  no  initial  angular  momentum.  Later, 
the  effects  of  hinging  at  a  support  will  be  explored.  As  an  example  of 
another  extreme,  initial  conditions  on  the  same  unit,  equal  to  a  speed  of 
5  ft/s  (1.52  m/s)  at  the  departure  time  of  54  ms  yielded  a  final  resting 
place  at  0.847  s  equal  to  a  downwind  displacement  of  10.4  ft  (3.17  m)  and  a 
crosswind  displacement  of  2.40  ft  (0.732  m) .  From  Fig.  14B,  such  a  result 
does  not  seem  to  have  been  likely  for  units  originally  located  in  the 
center  of  the  sidewall. 

Through  transport  calculations  of  this  kind  and  some  postshot  debris  dis¬ 
tribution  data,  it  appears  possible  to  deduce  the  sequence  of  wall  failure 
and  even  something  about  the  mode  of  failure. 

Plywood/Concrete  Cube  (DIAL  PACK) 

At  operation  DIAL  PACK, over  a  dozen  plywood  cubes  filled  with  concrete 
were  placed  on  the  ground  at  the  15,  30,  50,  and  100  psi  (103,  207,  345 
and  689  kPa)  contours,  and  final  downwind  displacements  recorded  (9). 

The  cubes  weighed  approximately  65  lbs  (29.6  kg)  and  were  1  ft  (.305  m) 
on  edge.  The  ground  was  not  smoothed  but  left  natural.  Judging  from 
their  final  attitudes,  the  cubes  travelled  by  tumbling  motion.  There  are 
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Trajectory  of  Center  of  Mass 


Figure  13A  Simulated  Trajectory  of  CMU  from  Middle  of  Sidewall  (Downwind) 
P.,,  =  30  psi  ( 207  k  Pa  ). 


no  motion  pictures  available,  so  far  as  is  known. 


Before  DEBRIS  could  be  used  in  the  simulation  of  tumbling  cubes,  it  was 
modified  to  account  for  attenuation  of  peak  blast  pressure  with  range 
from  ground  zero. 

In  the  simulation  at  15  and  30  psi  (103  and  207.  kPa)  the  cube  came  to 
rest  before  overturning  except  when  the  coefficient  of  friction  p  reached 
or  exceeded  1.0.  As  expected,  the  distance  travelled  declined  as  p  =  1.0, 
at  which  point  the  downwind  displacement  rose,  as  can  be  seen  in  Table  2. 


Table  2 


Coefficient 
of  Friction 


SIMULATED  PLYWOOD/ CONCRETE  CUBE 
Downwind  Displacement 

Displacement 


p 

15 

psi 

30 

psi 

(103. 

kPa) 

(207 

kPa) 

(ft) 

(m) 

(ft) 

(m) 

0.50 

15.3 

2.10 

0.75 

11.6 

3.54 

0.90 

10.1 

3.08 

41.0 

12.5 

0.95 

4.14 

1.26 

39.2 

11.9 

1.00 

>10.5 

3.20 

63.8 

19.4 

In  all  the  cases  presented  in  the  table,  airblast  pressure  drag  when  there 
was  no  tumbling  was  C^  =  1.05.  Ground  spring  constants  were 

K  .  =  2009  (3.5  x  108  )  and  K  =  .201  (3.5  x  104  ).  For 

gl  in  cm  g2  in  cm 

tumbling  cubes  or  cube-like  shapes  a  drag  coefficient  averaged  over  edge- 

on  and  head-on  attitudes,  i.e.,  C^  =  — 1 ^ 1 —  =  0.925,  was  used.  In 

every  case  the  drag  force  was  based  on  presented  area  and  base  drag  was 

considered  to  be  included  in  C, .  Within  the  limits  of  variation  of  K  . 

d  gl 

and  Kg2  explored,  the  influence  of  the  ground  spring  on  net  downwind  trans¬ 
lation  of  the  cube  appeared  small.  Definitive  work  on  the  effect  of  spring 
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constants  remains  to  be  done. 


The  trajectory  plotted  in  Figure  15  shows  clearly  why  the  net  downwind 
displacement  becomes  larger  for  p  =  1.0  than  for  p  <  1.0:  a  high  value 
of  friction  coefficient  can  produce  a  simulation  of  tumbling.  When  p  -  0.95 
the  trajectory  corresponding  to  that  of  Figure  15  is  flat. 

Figure  16,  summarizing  the  results  of  these  cube  observations  and  calcu¬ 
lations,  makes  it  clear  that  the  simple  coulomb  friction  model  is  inade¬ 
quate  in  cases  of  long  distance  translation.  For  studying  wall  and 
building  collapse,  however,  where  transport  is  over  relatively  short 
distances,  the  present  simulation  is  satisfactory.  In  Fig.  16  the  solid 
curve  is  drawn  through  the  average  of  the  values  of  net  downwind  displace¬ 
ment  of  cubes  observed  in  the  field  as  a  function  of  the  incident  over¬ 
pressure  at  their  original  location.  The  observed  values  themselves  are 
represented  by  solid  dots.  The  simulated  displacements  appear  as  crosses 
(  p  =0.95  )  and  crosses  in  circles  (  p  =  1.0  ).  Differences  between 
calculation  and  observation  increase  dramatically  with  incident  over¬ 
pressure,  suggesting  that  a  velocity-dependent  or  viscous  loss  model  may 
be  appropriate  (10).  In  fact.  Figure  16,  as  well  as  the  results  previous 
simulations,  indicate  that  in  the  short  range  p  =0.50  is  the  best  choice; 
viscous  forces  apparently  become  important  only  at  high  speeds. 


Trajectory  of  Center  of  Mass 


Figure  15  Calculated  Trajectory  of  Tumbling  Cube,  DIAL  PACK 
Pso-  15  psi  <  103  kPa  )  ,  yi  -  1.0. 


Figure  16  Simulated  Translation  of  Plywood  /  Concrete  Cube,  DIAL  PACK 
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IV  COLLAPSING  BUILDINGS 


The  source  of  debris  of  concern  to  civil  defense  is  the  partial  or  com¬ 
plete  collapse  of  buildings  after  an  airblast.  The  easiest  source  of  wall 
debris  to  analyze  with  tools  developed  to  this  point  is  probably  the 
masonry,  load-bearing  structure.  The  earliest  observation  of  the  collapse 
of  a  full-size,  four-wall  structure  of  this  kind  occurred  at  Operation 
TEAPOT  (11).  Partial  collapse  was  observed  at  5  psi  (34.5  kPa)  when  the 
sidewalls  apparently  fell  outward  due  to  poor  anchorage  against  outward 
motion  and  blast  entry  from  the  high  reflected  pressure  against  the 
front  of  the  building.  The  front  elevation  was  much  less  damaged.  Down¬ 
wind  translation  of  the  sidewall  debris  was  not  great,  indicating  relatively 
late  failure.  In  any  case,  this  early  observation  is  not  a  good  one  for 
present  purposes.  Of  more  interest  are  the  four-walled  structures  analyzed 
at  Operation  PRAIRIE  FLAT  (12).  These  were  of  two  kinds.  The  first  was 
an  approximately  cubical  structure  10  ft  (3.05  m)  on  edge  made  of  one  tier 
of  mortared  CMU  carrying  a  light  wood  roof  (Figure  10).  Although  some  of 
these  buildings  had  relatively  large  windows,  this  discussion  deals  only 
with  the  windowless  examples.  The  second  was  a  five-foot,  cubical  roofed 
structure  made  of  plastic  interlocking  blocks  glued  together  only  at  the 
corners  of  the  building  (Figure  17  and  18). 

CMU  Building,  Prairie  Flat 

Although  the  incipient  collapse  overpressure  of  the  CMU  Building  has  not 
been  calculated,  it  must  be  quite  low.  Exposed  at  two  locations,  i.e., 

9  and  30  psi  (62.0  and  207  kPa) ,  it  was  easily  and  completely  demolished  by 
the  blast  at  both  places  (13).  The  two  postshot  photographs  (Figures 
14A  and  14B)  are  the  only  data  found  so  far  that  describe  the  result  of 
these  experiments.  Any  postulated  sequence  of  collapse  events  based  on 
such  a  limited  amount  of  evidence  must  be  extremely  tentative. 
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Source:  Reference  12 


figure  17  Plastic  Block  Building,  PRAIRIE  FLAT. 
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PRAIRIE  FLAT  U/8  Scale) 


Source  : 


I - 1 - ^ — 

One  inch 


Reference  12 


Figure  18 


Plastic  Blocks,  PRAIRIE  FLAT. 


BRACOB,  short  for  Blast  Response  and  Collapse  of  Buildings,  is  a  computer 
code  designed  to  treat  simultaneously  blast  impact,  reflection  and  clearing, 
room-filling  or  backloading,  as  well  as  wall  deflection  and  collapse 
throughout  one  floor  of  a  building  (14).  As  the  simulated  blast  front 
sweeps  over  the  building,  the  appropriate  external  and  internal  wall  loads 
are  applied  as  functions  of  time.  Upon  collapse  of  a  wall,  the  interior 
pressures  against  certain  of  the  remaining  walls  are  adjusted  to  the  levels 
prevailing  outside  the  collapsed  wall.  At  this  stage  in  its  development, 
the  calculation  uses  an  arbitrary  delay  between  the  instant  of  wall  failure 
(i.e.,  when  it  can  no  longer  support  a  vertical  load)  and  the  time  of  trans¬ 
fer  of  the  exterior  pressure  to  the  interior.  The  magnitude  of  this  delay 
can  determine  the  final  disposition  of  wall  debris  and,  at  present,  there 
is  no  certain  way  of  calculating  it. 

Figure  19  presents  the  result  of  a  BRACOB  simulation  of  the  impact  of  the 
30-psi  (207-kPa)  blast  against  the  PRAIRIE  FLAT  CMU  building.  The  arbi¬ 
trary  delay  used  in  transferring  exterior  pressure  to  interior  was  4  ms.* 

In  the  figure  the  solid  line  represents  a  wall  before  deflection  and  the 
dashed  line  shows  schematically  the  deflection  of  the  wall  (either  inward 
or  outward).  The  deflection  is  deliberately  shown  as  zero  at  the  vertical 
wall  edges  and  greatest  at  the  center  to  represent  the  bowing  or  "arching" 
of  horizontally  loaded  walls  supported  at  the  edges.**  When  the  solid  line 
disappears,  the  wall  has  "collapsed",  i.e.,  the  time  delay  of  4  ms  after 
"failure"  has  passed. 

According  to  Figure  19,  collapse  of  the  front 


*  Values  of  other  structural  wall  parameters  during  this  simulation  were: 

•  Support  case  2 

•  Modulus  of  rupture,  12.9  psi  (88.9  kPa) 

•  In-plane  vertical  lead,  10.5  lb/in  (188  kg/m) 

•  Modulus  of  elasticity,  3  x  10^  psi  (2.07  x  10  ^  kPa) 

•  Density,  110.6  pcf  (1775  kg/nf*) 

•  Thickness,  7.63  (.194m). 

**  Photographic  data  exist  that  illustrate  the  cracking,  deflection,  and 

collapse  of  single  wall  panels  for  test  purposes.  Usually,  in  these  cases 


44 


luded  ) 
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wall  takes  place  at  12  ms  after  blast  arrival,  wti  •>.  m.  w.i  1  i 

are  still  moving  inward  and  the  rear  wall  is  h.ir.il-.  u  :  .•  :■■■..  .in.  . 

there  are  no  windows,  collapse  of  the  front  wal  1  pi.  vide-,  tin  1  it  i 
backloading  of  the  walls,  unless  blast  has  entered  through  t  • .»-  r..of  . 

At  19  ms  the  sidewalls  are  near  the  point  of  maximum  inw.ir..  d.  t  1.  .  t  ion. 
Under  normal  circumstances  they  would  now  begin  to  move  outwarn,  more 
or  less  reversing  (except  for  hysteretic  effects)  their  inward  path. 
However,  the  front  wall  is  now  "collapsed"  although  it  is  too  early  to 
see  the  opening  cleared  of  debris.  It  is  very  likely  11  ms  after  front 
wall  "failure"  there  is  a  substantial  amount  of  support  still  present  at 
the  leading  edges  of  the  sidewalls.  How  long  this  continues  is  not  known, 
but  the  subsequent  sequence  of  events  in  the  sidewalls  (as  calculated  by 
BRACOB  and  represented  in  Figure  19  for  times  greater  than  approximately 
40  ms)  does  not  seem  to  be  consistent  with  the  debris  pile  in  Figure  14B. 
After  44  ms  BRACOB  shows  the  sidewalls  arching  outward  as  if  fully  sup¬ 
ported  at  the  front  edges.  This  is  probably  unrealistic.  On  the  other 
hand,  BRACOB  also  shows  the  rear  wall  collapsing  at  39  ms.  If  this  does 
in  fact  happen  at  that  time,  the  blast  wind  (still  early  in  its  positive 
phase)  will  blow  through  the  structure  and  catch  the  inwardly  arched 
sidewalls . 

Simulations  with  DEBRIS  indicate  a  subsequent  translation  to  positions 
of  (a)  first  ground  contact  and  (b)  final  resting  place  in  the  general 


**  the  panel  interacts  in  some  way  with  a  rigid  or  nearly  rigid  frame  along 
its  edges.  One  such  test  was  performed  at  the  4.5-psi  (31.0-kPa)  contour 
at  Shot  Encore,  Operation  Upshot-Knothole.  Some  of  the  crack  patterns 
in  masonry  panels  seen  in  motion  pictures  from  Encore  are  sketched  in 
Figure  20  (15).  Fragments  are  varied  in  shape  but  there  are  generally  two 
classes:  those  that  are  "hinged"  along  an  edge  due  to  the  support  and 

those  that  are  punched  straight  out.  Usually,  the  latter  kind  are  found 
near  the  center  of  the  panel.  (In  examining  these  crack  patterns  it  must 
be  remembered  that  the  panels  tested  at  Encore  were  supported  in  an  en¬ 
tirely  different  manner  than  those  in  the  PRAIRIE  FLAT  CMU  house:  the 
failure  patterns  may  be  quite  different.  Note  also  that  the  arching 
drawn  in  Figure  19  represents  a  two-dimensional  effect,  i.e.,  the  figure 
cannot  show  the  bowing  that  may  be  occurring  in  a  vertical  plane.) 
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area  suggested  by  Figure  14B,  i.e.,  outward  and  downstream.  There  is 
some  dependence  on  the  exact  position  and  speed  of  the  inwardly  arching 
sidewall  fragment  at  the  moment  it  is  caught  by  the  through  flowing  blast 
and  on  the  degree  of  edge  support  at  the  front  wall  still  existing  at 
that  moment.  Assuming  all  support  is  lost  at  time  t  and  that  inward 
speed  and  deflection  are  and  ,  respectively,  then  displacements, 

AX  and  AY,  found  by  DEBRIS  are  shown  in  Table  3.  Displacements  to  first 
ground  contact  and  to  final  resting  place  are  reported.  These  results 
are  for  fragments  originally  receiving  suppo/t  at  the  vertical  edge 
(side  hinge)  on  the  GZ  side.  AY  and  AX  refer  to  crosswind  and  downwind 
travel,  respectively.  Z^  is  the  original  height  above  ground  of  the  CM. 
Predictions  by  DEBRIS  are  evidently  consistent  with  at  least  part  of  the 
debris  pattern  in  Figure  14B.  There  are,  of  course,  fragments  from  the 
portion  of  the  sidewall  that  hinges  at  the  vertical  edge  away  from  GZ. 

It  is  not  clear  when  this  material  is  exposed  to  the  flow-through  but  it 
is  reasonable  that  it  is  translated  to  the  rear  of  the  building  site. 

These  displacements  are  sketched  in  Figure  21. 

Disposal  of  the  front  wall  debris  is  also  important,  especially  if  it 

strikes  the  rear  wall  before  the  rear  wall  fails.  BRACOB  predicts  front 

wall  failure  at  8  ms  when  inward  speed  at  the  central  point  is  47  ft/s 

(14.3  m/s).  An  unhinged  fragment  not  in  ground  contact  originally  passes 

the  position  of  the  rear  wall  at  100  ms.  Since  this  is  the  first  kind 

of  front  wall  fragment  to  cross  the  depth  of  the  building,  the  interaction 

between  fragments  from  front  and  rear  walls  most  probably  takes  place 

well  downwind  of  the  building  site.  Incidentally,  first  ground  contact 

of  a  simulated  unhinged  front  wall  fragment  2  ft  x  2  ft  (.610  m.  x  .610  m) 

with  Z  =  5.0  ft  (1.52  m)  occurs  at  75  ft  (22.9  m) ,  and  at  319  ft  (97.2  m) 
c 

downwind  it  is  still  moving.  There  is  certainly  no  evidence  of  such  enor¬ 
mous  transport  in  Figure  14B.,  suggesting  the  importance  of  viscous  ground 
losses,  fragment-fragment  interaction,  or  errors  in  the  estimates  of 
departure  times  or  speeds. 
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Side-  and  bottom-hinged  fragments  from  the  front  wall  have  also  been 
simulated.  (There  appeared  to  be  essentially  no  structural  support  at 
the  roof,  thus  a  top-hinged  fragment  was  not  simulated).  A  simulated 
side-hinged  fragment  from  the  front  wall  originally  in  contact  with  the 
ground  is  flat  on  the  ground  at  a  downwind  displacement  equal  to  35.1  ft 
(10.7  m)  and  at  a  lateral  displacement  toward  the  hinge  equal  to  11.4  ft 
(3.47m).  While  still  partially  upright,  this  simulated  fragment  passes 
the  position  of  the  rear  wall  at  109  ms.  Its  lateral  displacement  at 
that  time  is  2.4  ft  (.731  m) .  At  that  same  time  a  simulated  bottom- 
hinged  fragment  from  the  rear  wall  is  less  than  5.2  ft  (1.58  m)  downstream 
and  a  simulated  side-hinged  fragment  near  the  rear  wall  is  displaced  less 
than  1  ft  (.304  m)  laterally  and  about  5  ft  (1.52  m)  downwind.  The  po¬ 
tential  for  debris  interaction  (and  greater  energy  dissipation  than  cal¬ 
culated  by  BRACOB)  is  thus  significant.  However,  it  would  appear  in  any 
case  that  the  initial  fragment  speeds  derived  from  BRACOB  are  excessive 
and  lead  to  unrealistically  large  downwind  transport  from  front  and  rear 
walls . 

The  data  in  the  9-psi  (62-kPa)  experiment  (Figure  14A)  are  less  clear  than 
those  taken  from  the  30-psi  (207-kPa)  experiment  (Figure  14B) .  For  one 
thing,  the  righthand  sidewall  in  Figure  14A  appears  partially  intact.  Also, 
evidently,  the  front  wall  debris  did  not  clear  the  building  site. 

A  BRACOB  simulation  of  the  PRAIRIE  FLAT  CMU  building  in  a  9-psi  (62-kPa) 
blast  predicts  wall  collapse  times  and  directions  as  follows*: 


Wall 

Time 

Direction 

front 

17.7  ms 

inward 

side 

30.0 

inward 

rear 

156.0 

outward 

At  the  time  the  side  walls  fail,  the  rear  wall  is  inwardly  deflected  1.5 
inches  (38  mm).  Thus,  if  the  simulation  is  at  all  close  to  reality,  there 
does  not  appear  any  likelihood  that  the  collapsed  side  walls  will  be  caught 
in  classical  blast  flowing  through  the  building.  Since  positive  overpressure 


*  Delay  for  debris  clearance  at  the  front  wall  was  taken  to  be  zero  in 
this  simulation 
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phase  duration  was  approximately  250  ms,  it  is  likely,  however,  that  the 
rear  wall  debris  went  downstream  of  the  building  and  such  a  disposition 
accords  with  Figure  14A.  Although  no  simulations  with  code  DEBRIS  of  this 
event  have  as  yet  been  made,  the  long  positive  phase  duration  probably  means 
that  enough  drag  exists  at  a  very  late  time  to  push  some  of  the  inwardly 
collapsed  side  wall  fragments  outward,  and  this  too  is  consistent  with  the 
observed  debris  distribution  on  the  left  side  of  the  building  in  Figure  14A. 

Plastic  Block  Model  Building 

These  models,  exposed  to  3  and  6  psi  (20.7  and  41.4  kPA)  at  PRAIRIE  FLAT, 
were  designed  to  provide  airblast  debris  similar  to  that  from  larger  masonry 
holdings.  No  BRACOB  simulations  of  collapse  have  been  attempted,  but  a 
cursory  examination  of  the  debris  data  reported  from  these  experiments  re¬ 
veals  a  pattern  of  distribution  similar  to  that  already  seen  in  the  CMU 
buildings.  A  tiny  sample  of  this  voluminous  data  is  reproduced  in  Figure  22, 
from  which  it  appears  clear  that  arching  in  these  models  took  place.  The 
glued  corners  created  realistic  support  conditions.  The  first  sketch  shows 
expected  lateral  displacement  due  to  vertical  edge  support.  The  kind  of 
rearward  and  lateral  displacement  in  the  second  sketch  seems  to  confirm  the 
expectation  of  early  blast  flow-through  while  the  sidewall  is  displaced 
inward  (as  postulated  in  Figure  19).  And  finally,  the  distribution  in  the 
third  sketch  suggests  an  outward  failure  of  the  rear  wall  while  vertical 
edge  support  from  the  sidewalls  still  existed. 
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Figure  22.  Observed  Debrjs  Distribution  from  Plastic  Block  Models 
PRAIRIE  FLAT,  Pso  -  3  psi  (  20.7  kPa  ) 
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V  ENHANCEMENT  OF  BRACOB 


The  present  effort  included  a  substantial  enhancement  of  BRACOB  in  pre¬ 
paration  for  the  calculation  of  urban  debris  distributions.  Storage 
was  increased  to  handle  large  buildings  and  capability  was  added  to 
treat  reinforced  concrete  and  masonry  walls  arching  in  a  yielding 
frame  (16). 

To  test  the  new  capabilities,  a  partial  response  to  a  30  psi-(207-kPa) 
head-on  blast  was  simulated  for  one  floor  of  the  Henry  R.  Landis  State 
Hospital  in  Philadelphia.  An  elevation  of  the  wing  that  was  studied 
is  shown  in  Figure  23.  A  typical  floor  plan  was  simplified  and  analyzed 
in  terms  of  an  orthogonal  grid  (Figure  24)  and  grid  intervals  taken 
from  construction  drawings  (Table  4).  Finally,  Figure  25  contains 
two  computer-drawn  floor  plans  of  the  typical  floor  of  the  hospital 
after  blast  arrival.  (Interior  doors  are  treated  as  closed.)  Collapse 
of  the  front  walls  at  this  pressure  is  almost  immediate. 
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HENRY  R.  LANDIS  STATE  HOSPITAL 


Table  IV 


Grid 


Dimensions,  Henry  R  Landis  State  Hospital. 


DIMX 

DIMY 

1  - 

1-829 

36- 

.609 

1  - 

5  080 

2  - 

.762 

37- 

1  .117 

2  - 

2  133 

3- 

1.117 

38- 

.254 

3  - 

203 

4  _ 

.254 

39- 

1  067 

4  - 

.203 

5- 

,508 

40- 

.609 

5  - 

2  146 

6- 

.559 

41  - 

.609 

6- 

3.580 

7  — 

.609 

42- 

1,067 

7- 

997 

8- 

.609 

43- 

254 

8- 

1  352 

9- 

1.067 

44- 

1.117 

9- 

.203 

10- 

.254 

45- 

609 

10- 

1  .829 

11  - 

1.117 

46- 

609 

11  - 

4  .140 

12- 

.609 

47- 

1.117 

12- 

1  .041 

13- 

.609 

48- 

.2  54 

14- 

1  .117 

49- 

1  067 

15- 

.254 

50- 

609 

16- 

.508 

51- 

609 

17- 

559 

52- 

1  067 

18  - 

.609 

53- 

.254 

19_ 

.609 

54- 

1.117 

20- 

1.067 

55- 

609 

21  _ 

.254 

56  - 

.609 

22- 

1.117 

57- 

1.117 

23  _ 

609 

58- 

.254 

24  - 

.609 

59- 

1  067 

25  _ 

.565 

60- 

.609 

26- 

.552 

61  - 

.609 

27  _ 

.254 

62- 

1.067 

28  _ 

1.067 

63- 

254 

29  _ 

609 

64- 

1.117 

30  - 

.609 

65- 

.609 

31  _ 

559 

66- 

609 

32  _ 

.508 

67- 

1.117 

33  - 

.254 

68- 

.406 

34  _ 

1.117 

69- 

.914 

35  _ 

.609 

70  - 

2.235 

ALL  DIMENSIONS  IN  METERS 
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VI  CONCLUSIONS  AND  RECOMMENDATIONS 

The  wall  collapse  (BRACOB)  and  debris  transport  (DEBRIS)  simulations 
appear  now  to  begin  serving  as  aids  in  understanding  the  fire-blast 
interaction.  Wall  failure  sequences  and  debris  locations  and  timing 
are  reasonably  well  simulated,  provided  transport  is  over  relatively 
short  distances.  However,  long-distance  transport,  even  of  heavy 
masonry  fragments,  does  not  appear  satisfactory  and  requires  investiga¬ 
tion  of  the  fragment  departure  velocities  and  a  possible  ground- 
fragment  interaction. 

In  addition  to  seeking  corrections  for  these  deficiencies,  it  is 
recommended  that  a  field  test  opportunity  be  sought  to  test  the 
theories  of  debris  creation  and  short-range  transport  from  wall 
masonry  buildings. 
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